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Abstract 

The pinning of flux lattices by weak impurity disorder is studied in the 
absence of free dislocations using both the gaussian variational method and, 
to 0(e = 4 — d), the functional renormalization group. We find universal 



logarithmic growth of displacements for 2 < d < 4: (u(x) — u(0)) 2 ~ log \x\ 

and persistence of algebraic quasi-long range translational order. When the 

two methods can be compared they agree within 10% on the value of A&. We 

compute the function describing the crossover between the "random manifold" 

regime and the logarithmic regime. This crossover should be observable in 

present decoration experiments. 
74.60.Ge, 05.20.-y 
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It has been argued on the basis of various elastic models for vortex lattices, such as 
Larkin's model of independent random forces acting on each vortex, that arbitrarily weak 
disorder destroys translational order below four dimensions ||. There is considerable dis- 
agreement, however, on the exact behaviour of the density-density correlation function 
Ck{ t ) — (Px(O)px(^)) even in the simpler case where dislocations are excluded. While 
direct extensions of Larkin's model predict exponential decay in d = 3 it has been 
pointed out that beyond the Larkin-Ovchinikov pinning length L c the lattice behaves col- 
lectively as an elastic manifold in a random potential with many metastable states, leading 
to a different power-law roughening of the lattice and to stretched exponential decay of 
Ck(j-) 0-0]. On the other hand, Flory-type arguments were proposed making explicit use 
of the periodicity of the lattice leading to logarithmic roughening ||. The role of disloca- 
tions at weak disorder above two dimensions is presently unsettled, but Bitter decoration 
experiments show remarkably large regions free of dislocations and provide a strong moti- 
vation for a better understanding of the elastic model. Other related pinned elastic systems 
such as charge density waves, magnetic bubbles, Wigner crystal, are under current active 
experimental study PH|JTT| ]. 

In this Letter we take into account both the existence of many metastable states and 
the periodicity of the lattice. We are primarily interested in the triangular Abrikosov lattice 
(d = 2 + 1). We also mention the case of d = 2 + (point vortices in thin films) or d = 1 + 1 
(lines in a plane). We show that in the absence of dislocations, the translational correlation 
function has a slow algebraic decay in dimension larger than two, and thus quasi-long range 
order persists. Two important length scales control the crossover towards this asymptotic 
decay, i) When the mean square of the relative displacement B(x) = |([w(x) — w(0)] 2 ) of two 
lines as a function of their separation x is shorter than the square of the Lindemann length 
l T = (u 2 ), the thermal wandering of the lines averages enough over the random potential and 
the model becomes equivalent to the random force Larkin model for which B(x) ~ |x| 4 ~ d . 
At low enough temperature, It is replaced by the superconducting coherence length £ (i-e. 
the correlation length of the random potential [fj-0]). ii) For l T < a 2 , B(x) ~ x 2v : 
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this is the random manifold regime where each line sees effectively an independent random 
potential, iii) For x > £, where £ corresponds to a relative displacement of the order of 
the lattice spacing a, B(x = £) ~ a 2 , the periodicity of the lattice becomes important. 
We find B(x) ~ Ajlog \x\ where Ad is a universal amplitude depending on dimension only 
and isotropy is recovered at large distances. This leads to quasi long range order CV (r) ~ 
(l/r) Ad . We have computed the full crossover function in d = 3 (Fig. 1). It suggests that all 
the above regimes could be observed by analysis of the dislocation-free decoration samples. 



These results are obtained using the Mezard-Parisi variational method [|12|] first applied in 
this context by Bouchaud, Mezard and Yedidia (BMY) ||[7|]- Our results are at variance 
with BMY, for reasons detailed below. In addition, we perform an e = 4 — d expansion using 
the functional renormalization group. The amplitudes obtained by these two rather 
different methods agree at order e within 10%. In d = 2, thermal fluctuations are important 
(It = oo) and the random manifold regime is much reduced. We find a modified Larkin 
regime with T-dependent exponents B (x) ~ |x| 2iy(T ) and a long distance logarithmic regime. 
Details can be found in [|T^]. 

We denote by Ri the equilibrium position of the lines labeled by an integer i, in the x — y 
plane, and by u(Ri, z) their in-plane displacements, z denotes the coordinate perpendicular 
to the planes. For weak disorder a/£ <C 1 it is legitimate to assume that u(Ri,z) is slowly 
varying on the scale of the lattice and to use a continuum elastic energy, as a function of the 
continuous variable u(x,z). Impurity disorder is modeled by a gaussian random potential 



V(x, z) with correlations: V(x, z)V(x', z') = A(x — x')5(z — z') where A(x) is a short range 
function of range £o an d Fourier transform A q . The total energy is: 

H e i = ^J d 2 xdz[(c u - c m )(d a u a ) 2 + c m (d a U/3) 2 + c 44 (d z u a ) 2 ] + J d 2 xdzV(x, z)p(x, z) (1) 

where a, (3 denote in-plane coordinates and the density is p(x, z) = J2i S(x — Ri — u(Ri, z)). 
Although we have also performed the calculations directly on the Hamiltonian ([I]) [|T^] it 
is more enlightening to use the following decomposition of the density that keeps track of 
the discreteness of the lines. In the absence of dislocations, generalizing ||14|| , one intro- 



- / d*xdzJ2[^d a u a a d p ul+Y: ^cos(K.(u a (x,z)-u\x,z)))] 



duces the slowly varying field <p(x,z) = x — u(<p(x, z), z). The density can be rewritten 
as p(x,z) = p det[d a (f)p] J2k e iK '^ x,z ^ ~ p (l - d a u a (<j>(x,z),z) + Ek^o e iK ' x p K (x)), where 

p K (x) = e -iK-u(<j>{x,z),z) ig 

the usual translational order parameter defined in terms of the 
reciprocal lattice vectors K, and po is the average density. 

Using the replica trick on (|I|) the disorder term gives — 1/(2T) J2 a ,b I d 2 xd 2 x'dzA(x — 
x')p a (x, z)p b (x', z). The above decomposition for the density leads to our starting model: 

Ha = / 10 £ G ->S( 9 , q,)u%q, q.) (2) 

0,6 Zi ir^O 

with G 1 = (c 4 4g^+c 66 g 2 )Pj / 3+(c4 4 g2+c 11 g 2 )Pf /3 in the case of (|T|), where P^ = 5 a(3 -q a q p /q 2 
and P^g = qalp/q 2 - To be rigorous (0) should be written in term of u((f>(x, z), z). This 
however has no effect on our results. It leads only to corrections of higher order in Vw which 
we neglect since we work here in the elastic limit a/^ < 1 [13|. For clarity we present the 
calculation for the isotropic version of (§) in d dimensions with Gq 1 = cq 2 5 a p where c is an 
elastic constant and p^Ax = A for all K (x = (x, z)). The results for flU) is presented at the 
end. A related single cosine model was studied by Villain and Fernandez [15] using a real 
space RG. Our results confirm and extend their analysis. 

One can get an idea of the effect of various terms in (|2]), by using arguments similar to 
T6|j . Assuming that u varies of ~ a over a length £, the density of kinetic energy is ~ c(a/£) 2 



. The long wavelength part of the disorder is Hj~ ~ / d xV(x)d a u a (x) ~ 

A l/2/£l+d/2 and 

Hq~K ~ / d d xV(x) cos(K • { x — u(x)) ~ /\}l 2 /^ d l 2 where K is the first reciprocal lattice 

/ , \l/(2-d) ( , \l/(4-d) 

vector. There are thus two length scales £ g ^o ~ a (era /Al and £ ~ a (c a /A) 

The g ~ component of the disorder is therefore relevant only for d < 2 and the second 
term in (||) can be dropped for d > 2. Higher Fourier components V qr ^K disorder the lattice 
below d = 4. 

We now look for the best trial Gaussian Hamiltonian Ho in replica space, of the form 
HI: 
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1 f d^Q 

Ho = 2J j2^y G ^^ Ua ^ ' ( 3 ) 

Defining the self-energy G~ b = cq 2 5 a b — &ab, and G~ 1 (q) = Z)bG~ 6 (g) we obtain by mini- 
mization of the variational free energy F var = F + (H — H ) jj the saddle point equations: 

GM) = — 2 , m = ^ kV ^ m) ^ 

cq z nT ^ 

B ab (x) = ±-([u a (x) - u b (0)] 2 ) (5) 
d d q 



T J —\ d {Gaa{q) + G bb (q) -2cos(qx)G ab (q)) 



(2tt) c 

where n is the number of components of u. 

The replica symmetric solution B a ^ b (x) = B(x), which mimics the distribution of dis- 
placements by a single Gaussian, is always unstable for 2 < d < 4. This is seen from the 
eigenvalue A of the replicon mode |12| : 



A = 1 - - E ^e"" 2T/ ^ Gc(P) / Al G2 M ( 6 ) 



Introducing a small regularizing mass in G c : G c (q)~ l = cq 2 + fi 2 we find, when \i — ► 0, 
that for d < 2 the replica symmetric solution is always stable and disorder is irrelevant. 
For d = 2 the condition becomes jU -2 ^ 1 4^ < 1 for small /i. Thus there is a transition at 
T = T c = Attc/Kq between a stable high-T phase where disorder is irrelevant and a low-T 
glassy phase where the symmetric saddle point is unstable. For 2 < d < 4 it is always 
unstable and disorder is always relevant. 

We now find a replica symmetry breaking solution for 2 < d < 4, the d = 2 case being 
discussed later. Following |T2[ we denote G(g) = G aa (q) and parametrize G ab (q) by G(q,v) 
where < t> < 1, and similarly for B ab (x). The saddle point equations become: 

ni K 

Expression (|]) and the algebraic rules for inversion of hierarchical matrices [12] give: 

, \ , N T c , f d d q 2Ta'(w) 
B(0,v) = B(0,v c ) + dw : . \ , „ 2 (8) 

./„ j (27r) d (Gc(g)- 1 + [(x](») 



where [o](v) = ua{y) — Jq dwa{w) and v c is the breakpoint such that a(v) is constant for 
v > v c . 5(0, v c ) is a nonuniversal quantity B(0, v c ) ~ £q + 

To discuss the large distance behaviour x ^> £ it is enough to keep K 2 = Kq in ([?[) 
since _B(0,t>) ^> a 2 . In that case, taking the derivative of (|7|) with respect to v, using 
[cr]'(f ) = va'{v) and (^) one finds the effective self energy: 

[a](v) = {v/vo) 2 ' 8 (9) 

and v = 2KlTc d c- d/2 /{4: - d). The energy fluctuation exponent is 9 = d — 2. Energy 
fluctuations are of order T/v and the large scale behaviour is controlled by small v. One 
can now compute the correlation functions: 

{(u(x) - u(0)r-) = 2nT J-0- d (l- cos(qx))G(q) (10) 



G(,) = i 5 ( 1+ / I *J^> (11) 
cq l Jo v z cq z + [a\ (v) q d 



with Z d = (4 - d)/{TK 2 S d ) and 1/S d = 2 d - 1 n d / 2 T[d/2}. Thus for 2 < d < 4 we find 
logarithmic growth: 



2n 



({u{x)-u{V)) 2 ) = — i A d \og\x\ (12) 
A o 

with A d = 4 - d. 

We now give the full crossover function for d = 3 for model (0) with PgA^ = A. 
The crossover length is f = a 4 c 2 /(27r 3 A). Defining /i(z) = £ P (A x /A^ )P 4 e- 2P , where 
K = 2irP/a, the solution is, in parametric form: 

v = v^h'(z) [a] 1 ' 2 = c^r'Kz) (13) 

with B(0,v) = a 2 z/(2ir 2 ) and = 27r 4 TA/(a 6 c 3 ) ~ It/C- The mean square displacement 
B is: 

B{x) = i 6(f ) &(*) = /o°° dt^^f{xh{t)) (14) 
/(a?) = 1-1(1 -e-) (15) 
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For x <C £, 5(0, v) is very small and h(z) ~ l/z n / 2+2 and therefore B{x) ~ x 2i/ with z/ = 1/6 
for n = 2. This corresponds to the random manifold regime W^\. The crossover function 



(0-0) was derived assuming <C v c , equivalent to It <C £• At scales such that B(x) is 
smaller than l T or £q one is in the regime v > v c and one recovers the replica symmetric 
propagator G(q) ~ 1/g 4 for q 2 ^> [cr](t> c ), and Larkin's model behavior. 

In the vortex lattice ([!]), shear deformations dominate (cqq <C Ch) and the crossover 
length becomes £ = c% 2 c l J^ a 4 / (2ti 3 A) . B a/3 (x) = B T P^(x) + B L P^ /3 (x) depends on the 
direction of x and we find: 

Bl*{x) = ^tAj) + #)] (16) 



where fr^r have an expression similar to (|14|) with / replaced by 



h{x) = \ ~ ~ 2 + (- + -^e"* / T (x) = /(*) - / L (s) (17) 
2 x z x x z 



The crossover function is shown in Fig. |1|. We find that if £/a ^> 1 all curves should 
scale when plotted in units of x/£. Such a crossover should be observable in decoration 
experiments. The ratio R = Bt{x)/Bl(x) crosses over at x = £ from R = 4/3 towards 
R = 1 as isotropy is restored at large scale. 

As d — > 2 + the function [cr](w) in @ vanishes for v < vq = TKq/(Attc). Thus in d = 2 
for T < T c there is a one-step replica symmetry broken solution with [<j](v) = for v < v c 



and [cr](i>) = [o"](i> c ) for v c < v < 1 |R|. In the glassy phase T < T c the variational 
method predicts the following two regimes. For T moderate, the short distance Larkin 
regime B(x) ~ x 2 crosses over directly towards the asymptotic logarithmic growth ( |T2|) at 
a length £ = a(<^a 2 /A) 1/(2_TK ° /27rc) . At low T, Tlog(£/a)/(ca 2 ) < 1, there is, in addition, 
an intermediate random manifold regime. In d = 1 + 1 the starting model (|2]) becomes 
exact due to the absence of dislocations and a RG calculation for T « T c [|T^,|Tj],|n|] finds 
-B(x) ~ log 2 (x) at large x. If one believes in the validity of this RG in a glassy phase, the 
simple gaussian ansatz does not give the exact long range behaviour, due to the importance 
of fluctuations in d = 2. However it gives T c exactly and captures correctly the crossover 



towards a slower logarithmic regime. Using RG we have shown [13| that in d = 2 the Larkin 



regime is in fact anomalous with a continuously variable exponent: 



TK 



B{x) ~ a 2 (x/0 2 ~^ (18) 

for x < £ ( in the low-T regime £ is replaced by another length). Model (|D will probably 
apply in d = 2 + at scales shorter than the distance between dislocations. It could also 
describe a polymerized membrane on a disordered substrate, with very high dislocation core 
energy. 

A previous application of the variational method by BMY ||[7] led to the conclusion which 
we believe is erroneous, that the fluctuations are enhanced at large distances. They applied 
the same method to a model in which each line sees a different disorder. This amounts 
to introduce an extra and unphysical disorder in the original model (|TJ) with correlations 
decaying as l/\Ri — Rj\ x - The long wavelength part of such a disorder dominates since large 
global translations of the lattice can improve the bulk energy, whereas in the physical model 
the gain in energy from the long wavelengths components can only come from surface terms 
and is thus irrelevant for d > 2, as shown in the discussion following (Q). Indeed for A — > 
the amplitude they obtain vanishes. In fact one can simplify the the saddle point equations 
of |||7|] by noting that the x dependence of B(x, u) in these equations is unimportant, up to 
higher order terms in Vu. The resulting local model (0) is simple enough to allow for the 
exact solution (|16|). 

To complement the variational method we perform a functional renormalization group 
calculation on the isotropic model. To simplify we take u to be a scalar field (n = 1) and 
set c = 1. One defines the replicated Hamiltonian as: 



-"imp 



±J2 J d d xA(u a (x) - u b (x)) (19) 

ab 



The function A(z) is periodic of period 1 (2tcKq = 1). The RG equations to order e = 4 — d 
have been derived by D.S. Fisher [HJ for the random manifold problem: 

^ = (e - 4C)A + C^A' + 1(A") 2 - A"A"(0) (20) 



A factor Sd has been absorbed into A. The periodicity of the function implies that the 
roughening exponent is ( = 0. This allows to obtain the fixed point function A* (z) in the 
interval [0,1]: 

A*(,) = ^(l-^(1-;,) 2 ) (21) 

Values for other z are obtained by periodicity. The fixed point is stable except for a constant 
shift. The linearized spectrum is discrete and the eigenvectors are Jacobi polynomials. This 
function has a nonanaliticity [|l9j A* ( - 4 - ) (0) = oo. To compute the correlation function T = TG 
one uses the RG flow equation: 

T{q, T, A) = e dl f (qe l , Te {2 ~ d)l , A(Z)) (22) 

choosing e r q = l/a allows to obtain perturbatively T(q) = — A"*(0)/g 4 = aS A ~ d 'e/ (36 Sdq d )- 
Thus we find, at order e, a logarithmic growth ([12]) of line displacement with A^rg = 
e(27r) 2 /36 = l.lOe. This compares within 10% with A^var — £ which slightly underestimate 
fluctuations. For 2 < d < 4 the real space RG method of |I5| also predicts a log but does 



not allow to compute the universal prefactor A^ or the crossover function. The agreement 
between these methods, none being rigorous, lends credibility to the additional results in 
d = 3 obtained using the variational method. 

To conclude, we have shown that due to the periodicity of the lattice, the pinning by 
impurities becomes less effective at large length scales. This has consequences on the trans- 
port properties. In the flux creep regime energy barrier arguments @, and the exponents 
found here, lead to a non linear voltage-current relation V ~ exp[— l/(Tj M )] where fi crosses 
over from \i w 0.7 — 0.8 to /x = 1/2 as j decreases. 

We thank D.S. Fisher, T. Hwa and D. Huse for interesting discussions. Both authors 
thank AT&T Bell Laboratories, where part of this work was completed, for support and 
hospitality. PLD acknowledges support from NSF Grant 4-20508. 
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FIGURES 

FIG. 1. Plot of bi versus log|x/£| ( solid line ), as defined in the text, for the Abrikosov 
triangular lattice. The dashed line is the random manifold result. 
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